args Controls

********************************************************************************
* Wage distribution (UQR)
********************************************************************************

replace Creativity = creativityt_psd 


* 25th
foreach a of numlist 23 33 42 46 50 55{
qui rifhdreg lwage_I`a' `Controls' [pweight=1/pObs_I`a'], rob rif(q(25))
est sto e`a'
}

esttab e* using results/Tables/R_OLS_winter.tex, ///
replace se star(* 0.1 ** 0.05 *** 0.01) b(%9.3f) se(%9.3f) nonumbers collabels(none) ///
keep(Creativity CogAbility) ///
refcat(Creativity "25th percentile of log wages" " ", nolabel) ///
order(Creativity) ///
varlabels(Creativity "Creativity" CogAbility "Cognitive ability") ///
cells(b(star fmt(3)) se(par) ".") nogaps   alignment(c) frag ///
mtitles("Age 23" "Age 33" "Age 42" "Age 46" "Age 50" "Age 55") ///
posthead(\cline{2-7})  noobs   
est drop _all


* 50th
foreach a of numlist 23 33 42 46 50 55{
qui rifhdreg lwage_I`a' `Controls' [pweight=1/pObs_I`a'], rob rif(q(50))
est sto e`a'
}


esttab e* using results/Tables/R_OLS_winter.tex, ///
append `text' se star(* 0.1 ** 0.05 *** 0.01) b(%9.3f) se(%9.3f) nonumbers collabels(none) ///
keep(Creativity CogAbility) ///
refcat(Creativity "50th percentile of log wages" " ", nolabel) ///
order(Creativity) ///
varlabels(Creativity "Creativity" CogAbility "Cognitive ability") ///
cells(b(star fmt(3)) se(par) ".") nogaps   alignment(c) ///
posthead("") nomtitles frag   noobs  
est drop _all

* 75th
foreach a of numlist 23 33 42 46 50 55{
qui rifhdreg lwage_I`a' `Controls' [pweight=1/pObs_I`a'], rob rif(q(75))
est sto e`a'
}

esttab e* using results/Tables/R_OLS_winter.tex, ///
append `text' se star(* 0.1 ** 0.05 *** 0.01) b(%9.3f) se(%9.3f) nonumbers collabels(none) ///
keep(Creativity CogAbility) ///
refcat(Creativity "75th percentile of log wages" " ", nolabel) ///
order(Creativity) ///
varlabels(Creativity "Creativity" CogAbility "Cognitive ability") ///
cells(b(star fmt(3)) se(par) ".") nogaps   alignment(c) ///
posthead("") nomtitles frag  ///
stats(N, fmt(%11.0gc) labels("Observations"))
est drop _all
